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We study the formation of black holes by spherical domain wall collapse as seen by an asymptotic 
observer, using the functional Schrodinger formalism. To explore what signals such observers will 
see, we study radiation of a scalar quantum field in the collapsing domain wall background. The 
total energy flux radiated diverges when backreaction of the radiation on the collapsing wall is 
ignored, and the domain wall is seen by the asymptotic observer to evaporate by non-thermal "pre- 
Hawking radiation" during the collapse process. Evaporation by pre-Hawking radiation implies that 
an asymptotic observer can never lose objects down a black hole. Together with the non-thermal 
nature of the radiation, this may resolve the black hole information loss problem. 



I. INTRODUCTION 

Black holes embody the long-standing theoretical chal- 
lenge of combining general relativity and quantum me- 
chanics, with various proposals being advocated over the 
years to resolve paradoxes associated with black hole for- 
mation, evaporation and information loss. Resolution of 
these issues has become even more timely with the pos- 
sible formation and evaporation of black holes in particle 
accelerators in the framework of higher dimensional mod- 
els that have recently garnered much attention. The pro- 
cess of black hole formation is generally discussed from 
the viewpoint of an infalling observer. However, in all 
physical settings it is the viewpoint of the asymptotic 
observer that is relevant. More concretely, if a black hole 
is formed in the Large Hadron Collider, it has to be ob- 
served by physicists sitting on the CERN campus. The 
physicists have clocks in their offices and they watch the 
process of formation and evaporation in this coordinate 
frame. They must address questions such as: At what 
time did a black hole form? Is any information lost into 
the black hole? How long did it take for the black hole to 
evaporate? What is the spectrum of the decay products? 

The process of gravitational collapse has been studied 
extensively over the last few decades, from many differ- 
ent viewpoints, including 1-1-1 dimensional models and 
modifications of general relativity [e.g. see Q). Unlike 
a large subset of this work, our analysis is in 3-1-1 di- 
mensions and within conventional general relativity. We 
model the general problem by choosing to study a collaps- 
ing spherical shell of matter, more specifically a vacuum 
domain wall. The physical setup of the problem and the 
fimctional Schrodinger formalism are described in Sec.HIl 

A crucial aspect of our analysis is that we address 
the question of black hole formation and evaporation 
as seen by an asymptotic observer. Initially, when the 
domain wall is large, the spacetime is described by the 
Schwarzschild metric, just as for a static star. From here 
on, the wall and the metric are evolved forward in time. 



always using the Schwarzschild time coordinate. We em- 
phasize that all our discussion, unless explicitly stated, 
refers to the Schwarzschild time, t, and this defines the 
time slicing of the spacetime. As is well known, the 
Schwarzschild coordinate system breaks down at a black 
hole horizon, and there is danger that our analysis will 
also break down at some point during the gravitational 
collapse. However, we do not encounter any such difficul- 
ties, suggesting that our calculation is self-consistent. A 
second danger is that the coordinate system may provide 
an incomplete description of the gravitational collapse 
spacetime. This remains a possibility. However, wc find 
that Schwarzschild coordinates are sufficient to answer 
the very specific set of questions we ask from the asymp- 
totic observer's viewpoint. Namely, does the asymptotic 
observer see objects disappear into a black hole in the 
time that he sees the collapsing body evaporate? And, is 
the spectrum of the radiation received ever truly thermal 
(even in the semiclassical approximation)? 

In Sec. lIIII we verify the standard result that the forma- 
tion of an event horizon takes an infinite (Schwarzschild) 
time if we consider classical collapse. This is not 
surprising and is often viewed as a limitation of the 
Schwarzschild coordinate system. To see if this result 
changes when quantum effects are taken into account, we 
address the problem of quantum collapse using a minisu- 
perspace version of the functional Schrodinger equation 
[31 in Sec. |TVl We find that even in this case the black 
hole takes an infinite time to form, contrary to some spec- 
ulations in the literature Q . 

In Sec.|V]wc consider the possible radiation associated 
with the collapsing shell by considering the interaction 
of a quantum scalar field and the classical background 
of a collapsing domain wall. We treat the problem us- 
ing the functional Schrodinger picture, which we relate to 
the standard Bogolubov treatment carried out in Sec. lVIl 
Here we find that the shell, even as it collapses, radiates 
away its energy in a finite amount of time. With some 
assumptions about the metric close to the incipient hori- 
zon, we conclude that the evaporation time is shorter 
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than what would be taken by objects to fall through a 
black hole horizon. This leads us to the conclusion that 
the asymptotic observer will sec the evaporation of the 
collapsing shell before he can see any objects disappear. 

We discuss our results from the point of view of an in- 
falling observer in Sec. IVIH where we attempt to recon- 
cile the fact that such an observer will not see substantial 
radiation with the observations made by an asymptotic 
observer. Our conclusions are summarized in Sec. IVIIIl 
where we elucidate a possible resolution of the informa- 
tion loss problem suggested by our results, together with 
a discussion of possible loopholes and future directions. 



II. SETUP AND FORMALISM 

To study a concrete realization of black hole formation 
we consider a spherical Nambu-Goto domain wall that is 
collapsing. To include the possibility of (spherically sym- 
metric) radiation we consider a massless scalar field, $, 
that is coupled to the gravitational field but not directly 
to the domain wall. The action for the system is 



S = dt^Xyf—g 
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(1) 



where the first term is the Einstein-Hilbert action for the 
gravitational field, the second is the scalar field action, 
the third is the domain wall action in terms of the wall 
world volume coordinates, (a = 0, 1,2), the wall ten- 
sion, a, and the induced world volume metric 



lab = g^^.daX'^dbX'' 



(2) 



The coordinates X^^ (^° ) describe the location of the wall 
and Roman indices go over internal domain wall world- 
volume coordinates C", while Greek indices go over space- 
time coordinates. The term iSobs m Eq. ^ denotes the 
action for the observer. 

We will begin first with the Wheeler-de Witt equation 
in order to explore and contrast quantum vs classical col- 
lapse of the domain wall, but we will eventually utilize 
the functional Schrodingcr formalism to study both col- 
lapse and radiation in this system. 

The Wheeler-de Witt equation for a closed universe is 



= 



(3) 



where H is the Hamiltonian and g^^, $, O] is the 

wave-functional for all the ingredients of the system, in- 
cluding the observer's degrees of freedom denoted by O. 
Note that the wave-functional, 4", is a functional of the 
fields but not of the spacetime coordinates. We will sep- 
arate the Hamiltonian into two parts, one for the system 
and the other for the observer 



(4) 



Any (weak) interaction terms between the observer and 
the wall- metric-scalar system are included in ifgys- The 
observer is assumed not to significantly affect the evolu- 
tion of the system and similarly for the system vis a vis 
the observer. The total wave-functional can be written 
as a sum over eigenstates 



;sys,t)vl/^bs(^?,*) 



(5) 



where k labels the eigenstates, cu are complex coeffi- 
cients, and we have introduced the observer time, t, via 



.9* 



obs 



dt 



^^obs*obs 



(6) 



With the help of an integration by parts, and the fact 
that the total wave-functional is independent of the 
Wheeler-de Witt equation implies the Schrodinger equa- 
tion 



' dt 



(7) 



For convenience, from now on we will denote the system 
wave-function simply by 'J and drop the superscript k 
and the subscript "sys" . Similarly H will now denote 
ffsys, and the Schrodinger equation reads 



(8) 



A general treatment of the full Wheeler-de Witt equa- 
tion is very difficult and we shall utilize the frequently 
employed strategy of truncating the field degrees of free- 
dom to a finite subset. In other words, we will consider a 
minisuperspace version of the Wheeler-de Witt equation. 
As long as we keep all the degrees of freedom that are 
of interest to us, this is a useful truncation. With this 
in mind, we only consider spherical domain walls and as- 
sume spherical symmetry for all the fields. So the wall is 
described by only the radial degree of freedom, R{t). The 
metric is taken to be the solution of Einstein equations 
for a spherical domain wall. The metric is Schwarzschild 
outside the wall, as follows from spherical symmetry Q 

ds^ = -{l-^)dt^ + {l-^)-^dr^ +r^dn'' , r > i?(t) 
r r 

(9) 

where, Rs ~ 2GM is the Schwarzschild radius in terms 
of the mass, M, of the wall, and 

dn^ = de^ + sin^ ^-^-^^ 



(10) 

In the interior of the spherical domain wall, the line ele- 
ment is flat, as expected by Birkhoff's theorem, 

r < R{t) (11) 



-dT^ + dr'^Wde'^W sir 2 o^^"^ 



ds^ 



The interior time coordinate, T, is related to the observer 
time coordinate, i, via the proper time, r, of the domain 
wall. 
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dR 
d7 



(12) 



and 



dt 



1 

B 



dR 
d^ 
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where 



R 



B = 1 

The ratio of these equations gives 



dT _ {l + B:if/^B 
It ^ {B + i?2)l/2 



B 



B 



1/2 



(13) 



(14) 



(15) 



where Rr = dR/dr and R = dR/dt. Integrating Eq. ([T5]) 
stiU requires knowing R{t) (or R{t)) function of r 
(ort). 

Since we are restricting our minisuperspace to fields 
with spherically symmetry, we need not include any 
other metric degrees of freedom. The scalar field can 
also be truncated to the spherically symmetric modes 
($ = $(t,r)).^ 

By integrating the equations of motion for the spherical 
domain wall, Ipser and Sikivie [3] found that the mass is 
a constant of motion and is given by 



1, 



M = -[v/l + i?? + + RV\ATTaR^ 



(16) 



where it is assumed that max(i?) > l/47rGcr Ql- This 
expression for M is implicit since Rs = 2GM occurs in 
B. Solving for M explicitly in terms of Rr gives 



M = A-KcrR^yi + i?2 _ 2'KGaR] 



(17) 



and with the relations between T and t given above we 
get 



M = inaR-" 



1 



2'KGaR 



(18) 



where Rt denotes dR/dT. 

We now discuss the classical collapse of the domain 
wall. 



III. CLASSICAL TREATMENT OF DOMAIN 
WALL COLLAPSE 

A naive approach to obtaining the dynamics for the 
spherical domain wall is to insert the spherical ansatz for 
both the wall and the metric in the original action. This 
would lead to an effective action for the radial coordinate 
R{t). However, it is known that this approach does not 
give the correct dynamics for gravitating systems. We 
find that this approach does not straightforwardly lead 
to mass conservation as given in Eq. (|16p . So we take 
the alternative approach of finding an action that leads 



to the correct mass conservation law. The form of the 
action can be deduced from Eq. ^TE\\ quite easily 



A-Kcj / dTR^ 



I- R^- 27rGaR 



(19) 



which can now be written in terms of the external time t 
Seff = -47rCT / dtR'^ 



^eff 

- 2nGaR^B 
and the effective Lagrangian is 



'B 



i?2 



B 



(20) 



-AnaR-^ 



'B 



R? 
B 



2TrGaR\ B 



(21) 

The generalized momentum, H, can be derived from 
Eq. (Ell) 



n 



AnaR^R 1 2TrGaR{l - B) 

VB [^B^-R2 ^B^-{1-B)R\ 

The Hamiltonian (in terms of R) is 
H = AiraB^I^R^ 



(22) 



2TrGaR 



Vb^-R-^ Jb^ - {l-B)m 



(23) 

To obtain iJ as a function of (i?,n), we need to elim- 
inate R in favor of 11 using Eq. (EH). This can be done 
but is messy, requiring solutions of a quartic polynomial. 
Instead we consider the limit when R is close to Rs and 
hence B — *■ 0. In this limit the denominators of the two 
terms in Eqs. (also in ([^5]) ) are equal and 



n 



47r/ii?2i? 



where 



and 



H 



^/bVb^ - ie 

^ = cr(l - 2nGaRs) 
AiT^iB^/^R^ 



(24) 



(25) 



(26) 



- 7?2 

= [{Bllf +B{4TTfiRY]^^^ (27) 

The Hamiltonian has the form of the energy of a rela- 
tivistic particle, ^/p^ + m? , with a position dependent 
mass. 

The Hamiltonian is a conserved quantity and so, from 

Eq. jig). 



B3/2j?2 



= h 



(28) 
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where h = H/Attij, is a constant. (Up to the approxima- 
tion used to obtain the simpler form of the Hamiltonian 
in Eq. ([26|) . the Hamiltonian is the mass.) 
Solving Eq. for R wc get 



R = ±B 1- 



1/2 



which, near the horizon, takes the form 
R 



1 BR^ 
±5 1-- 



2 /i2 



(29) 



(30) 



since i? — > as i? ^ Rs- 

The dynamics for R ^ Rs can be obtained by solving 
the equation R = ±B. To leading order in i? — Rs, the 
solution is 



R{t) « i?s + {Ro - Rs)e 



±t/Rs 



(31) 



where Rq is the radius of the shell at i = 0. As we 
are interested in the collapsing solution, we choose the 
negative sign in the exponent. This solution implies that, 
from the classical point of view, the asymptotic observer 
never sees the formation of the horizon of the black hole, 
since R{t) = Rs only as t ^ oo. This result is similar to 
the well-known result (for example, see [H]) that it takes 
an infinite time for objects to fall into a pre-existing black 
hole as viewed by an asymptotic observer . In our case 
there is no pre-existing horizon, which is itself taking an 
infinite amount of time to form during collapse. To see 
if this conclusion will change when quantum effects are 
taken into account {e.g. Sec. 10.1.5 of 0]) we now explore 
the quantum dynamics of the collapsing domain wall. 



IV. QUANTUM TREATMENT OF DOMAIN 
WALL COLLAPSE 

The classical Hamiltonian in Eq. ([?7|) has a square root 
and so we first consider the squared Hamiltonian 



^ BH Bn + B{ATT^R^ f 



(32) 



where wc have made a choice for ordering B and H in 
the first term. In general, we should add terms that 
depend on the commutator [B, H]. However, in the limit 
R — > Rs, we find 

Estimating H by the mass, M , of the domain wall, the 
terms due to the operating order ambiguity will be neg- 
ligible provided 



1 

'r's 



M 



where mp is the Planck mass. Hence the operator order- 
ing ambiguity can be ignored for domain walls that are 
much more massive than the Planck mass. 



Now we apply the standard quantization procedure. 
We substitute 



H 



.d_ 

'dR 



in the squared Schrodinger equation, 



Then 



To solve this equation, define 



u = i? -f- i?s In 



which gives 



SH 



R_ 
'r's 



du 



Eq. ([34]) then gives 



92^ 92* 



B{Att^R'^)H = 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



This is just the massive wave equation in a Minkowski 
background with a mass that depends on the position. 
Note that R needs to be written in terms of the coordi- 
nate u and this can be done (in principle) by inverting 
Eq. ([5^ . However, care needs to be taken to choose the 
correct branch since the region R G (Rs, 00) maps onto 
u G (—00, -l-oo) and R G (0, Rs) onto u G (0, —00). 

We are interested in the situation of a collapsing wall. 
In the region R ^ Rs, the logarithm in Eq. ([5^ domi- 
nates and 

R^ Rs + i?se"/-"^ 

We look for wave-packet solutions propagating toward 
Rs, that is, toward u — > — 00. In this limit 



B 







and the last term in Eq. p8|) can be ignored. Wave 
packet dynamics in this region is simply given by the free 
wave equation and any function of light-cone coordinates 
(u ± t) is a solution. In particular, we can write a Gaus- 
sian wave packet solution that is propagating toward the 
Schwarzschild radius 



1 



-(ti+t)^/2s^ 



27rs 



(39) 



where s is some chosen width of the wave packet in the 
u coordinate. The width of the Gaussian wave packet 
remains fixed in the u coordinate while it shrinks in the R 
coordinate via the relation dR = Bdu which follows from 



5 



Eq. ((36)) . This fact is of great importance, since if the 
wave packet remained of constant size in R coordinates, 
it might cross the event horizon in finite time. 

The wave packet travels at the speed of hght in the u 
coordinate - as can be seen directly from the wave equa- 
tion Eq. ([38]) or from the solution, Eq. ([39| . However, 
to get to the horizon, it must travel out to u = —oo, 
and this takes an infinite time. So we conclude that the 
quantum domain wall does not collapse to Rs in a finite 
time, as far as the asymptotic observer is concerned, so 
that quantum effects do not alter the classical result that 
an asymptotic observer does not observe the formation 
of an event horizon. 

The above analysis shows that the collapsing wall takes 
an infinite time to reach R = Rs- The analysis leaves 
room for processes by which the wave packet can jump 
from the {Rs, oo) region to the {0,Rs) region, without 
ever going through Rs- Note that this is different from 
tunneling through a barrier. In that case, the wave func- 
tion is non-zero within the barrier, and a small part of it 
leaks over to the other side of the barrier. In the present 
case, Rs occurs at w = — oo and so, if there is any bar- 
rier, it is infinitely far away. If there is to be a jump from 
outside to inside Rs-, it does not show up in the present 
description using the Wheeler-dc Witt equation. 

We have obtained the massive wave equation, Eq. 
by first squaring the classical Hamiltonian, Eq. (P7)) . This 
procedure eliminated the square root occurring in the 
Hamiltonian. It is possible that some other quantization 
procedure will yield different conclusions. In this context, 
we note, in fact, that we need not square the Hamiltonian 
to get rid of the square root provided we work in the near 
horizon limit. In that case 

H = [{Bnf + BiAnfiR^)'^] ^'"^ w ±BH (40) 

where the sign is chosen to make B non-negative. Then 
the Schrodinger equation again yields wave packets prop- 
agating at the speed of light in the [t, u) coordinate sys- 
tem and with the horizon located at it = — oo. 



V. RADIATION - SEMICLASSICAL 
TREATMENT 

If an external observer never sees the formation of an 
event horizon, we need to explore what radiation might 
be observed that characterizes gravitational collapse. To 
do so we consider a quantum scalar field in the back- 
ground of the collapsing domain wall. We do not con- 
sider gravitational radiation since this is excluded by our 
restriction to spherically symmetric modes in minisuper- 
space. In this section, we approach the problem using the 
functional Schrodinger equation since (i) we have already 
set up this approach and used it in the previous section, 
(ii) we believe the approach is more suited to treating the 
backreaction problem, and (iii) it allows us to calculate 
the total radiation of which Hawking radiation may only 
be a subset. To connect with earlier work, we discuss 



the problem of Hawking radiation using the conventional 
Bogolubov transformations in Sec. IVII 
The action for the scalar field is 



(41) 



We decompose the (spherically symmetric) scalar field 
into a complete set of real basis functions denoted by 



<i> = ^afc(i)/fe(r) 



(42) 



The exact form of the functions /fe (r) will not be impor- 
tant for us. We will be interested in the wavefunction for 
the mode coefficients {afc}. 

In the functional Schrodinger picture, we wish to 
find the wavefunctional by solving a functional 

Schrodinger equation. This is equivalent to finding the 
wavefunction of an infinite set of variables, V'({ofc},i), 
by solving an ordinary Schrodinger equation in an infi- 
nite dimensional space. The procedure (detailed below) 
is to find independent eigenmodes, for the system 

for which the Hamiltonian is a sum of terms, one for 
each independent eigenmode. Then the wavefunction fac- 
torizcs and can be found by solving a time-dependent 
Schrodinger equation of just one variable. 

Since the metric inside and outside the shell have dif- 
ferent forms, we split the action into two parts 



where 



Sin ^2Tr I dt J dr 



t 



2tt I dt I dr 

lR(t) 



1 - Rs/r 



(43) 



(44) 



1 - ^ ) {dr^f 



T is given by Eq. p^ . which with Eq. (P^)) . gives 



■ dT 

T = — ^ B 
dt 



n4' 



1/2 



(45) 



(46) 



As i? — > Rs, T ^ B ^ 0. Therefore the kinetic term in 
5in diverges as (i? — Rs)~^ in this limit and dominates 
over the softer logarithmically divergent contribution to 
the kinetic term from 5out • Similarly the gradient term in 
Sin vanishes in this limit and is sub-dominant compared 
to the contribution coming from S'out- Hence, 



S 



27r / dt 



1 



-R, 

BJo 
dr r^ I 1 



Rs 



dr r^{dt<i>)^ 
Rs 



(47) 
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where we have changed the hmits of the integrations to 
Rs since we are working in the regime R{t) ~ i?5. This 
approximation is vahd provided the contribution from 
r £ {Rs, R(t)) to the integrals remains subdominant, and 
also the time variation introduced by the true integration 
limit {R{t)) can be ignored. These requirements are not 
arduous. 

Now, we use the expansion in modes in Eq. (|42p to 
write 



S= dt 



(48) 



where M and N are matrices that are independent of 
R(t) and are given by 

i-Rs 

Mfe,.'=4^/ dr r^fk{r)fk'{r) (49) 



Nkk' = 47r / dr r 
iRs 



"l~ — \mf'k'{r) (50) 



Using the standard quantization procedure, the wave 
function ipiak,t) satisfies 



1 



where 



R 



Uk{M-^)kk'nk' + lakNkk'ak 



^ = i 



dt 
(51) 



Hfc = -I 



dak 



(52) 



is the momentum operator conjugate to aj,. 

So the problem of radiation from the collapsing domain 
wall is equivalent to the problem of an infinite set of 
coupled harmonic oscillators whose masses go to infinity 
with time. Since the matrices M and N are symmetric 
and real [i.e. Hermitian), it is possible to do a principal 
axis transformation to simultaneously diagonalize M and 
N (see Sec. 6-2 of Ref. [2] for example). Then for a single 
eigenmode, the Schrodingcr equation takes the form 



1 



Rs_ 
R 



2m 962 



1 



K¥ 



(53) 



where m and K denote eigenvalues of M and N, and b 
is the eigenmode. 

We re-write Eq. ([53]) in the standard form 



1 92 m 2/ N,2 



d 



where 



and 



K 



dt\l- 



1 



R^ 
R 



ml-Rs/R l-Rs/R 



(54) 



(55) 



(56) 



1 - Rs/R 



FIG. 1: Model for R{t). 

We have chosen to set ry(f = 0) = 0. 

To proceed further, we need to choose the background 
spacetime i.e. the behavior of R{t). The classical solu- 
tion in Eq. tells us that 1-Rs/R^ exp{-t/Rs) at 
late times. We are mostly interested in the particle pro- 
duction during this period. At early times, the behavior 
depends on how the spherical domain wall was created 
and we are free to choose a behavior for R{t) that is con- 
venient for calculations and interpretation. To be able 
to interpret particle production at very late times it is 
easiest to have a static situation. This can be obtained if 
we artificially take the collapse to stop at some time, tf. 
Eventually we can take t j ^ 00 to go over to the eternal 
collapse case. So our choice for R will be 



El 
R 



1 



(,-tf/Rs 



t e (-00,0) 
t e (i/, 00) 



(57) 



This choice does have the issue that the derivative of R 
has discontinuities at t = and t = tf. However, we 
shall show below that these discontinuities do not affect 
particle production. 

With the chosen behavior of R, the spacetime is static 
at early times and the initial vacuum state for the modes 
is the simple harmonic oscillator ground state, 

, ,1/4 

Hb,V = 0) = 



(58) 



at later times is 



Then the exact solution to Eq 

/ \ 1/4 
V;(6,,y) =e*"(''M -^j cxp 

. . ^(^^^ 

where p,, denotes derivative of p{rj) with respect to ry, and 

p is given by the real solution of the ordinary (though 

non-linear) differential equation 



Pr^,^+Uj'^{l])p = 4t 



with initial conditions 
P(0) = 



PviO) = 



(60) 



(61) 
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The phase a is given by 



air-j) 



_drf_ 
7W) 



(62) 



In Appendix [A] we discuss the behavior of p as — > Rs 
it — > oo). Also note that the solution for p and is 
continuous. 

Consider an observer with detectors that are designed 
to register particles of different frequencies for the free 
field <f> at early times. Such an observer will interpret the 
wavefunction of a given mode b at late times in terms 
of simple harmonic oscillator states, {^Pn}, at the final 
frequency, 



ujQe 



tf/2B.s 



(63) 



The number of quanta in eigenmode b can be evaluated 
by decomposing the wavefunction (Eq. I|59p ) in terms of 
the states, {</5n}, and by evaluating the occupation num- 
ber of that mode. To implement this evaluation, we start 
by writing the wavefunction for a given mode at time 
t > tf in terms of the simple harmonic oscillator basis at 
t = 0. 



where 



db ^im{b,t) 



(64) 



(65) 



which is an overlap of a Gaussian with the simple har- 
monic oscillator basis functions. The occupation number 
at eigenfrequency uj {i.e. in the eigenmode b) by the time 
t > tf., is given by the expectation value 



(66) 



In Appendix |B] we evaluate the occupation number in 
the eigenmode b and the result is given in Eq. (|B12p 



Nit, Co) 



1 - 



U)p 



(67) 



for t > tf. 

By calculating N, it can be checked that N remains 
constant for t < and also t > tf. Hence all the particle 
production occurs for < t < tf . There is a possibility 
that the particle production is due to discontinuities in 
the derivative of R at t ~ 0,t/. However, as we shall 
see below, the particle number grows with increasing tf, 
while the discontinuity at t = is fixed, and that at t = 
tf gets weaker. This indicates that particle production 
occurs only during < t < tf and is a consequence of 
the gravitational collapse. 

Now we can take the tf ^ oo limit. In Appendix \X\ we 
have shown that p remains finite but —oo as t > 



tf oo, provided loq 0. However, we are interested 
in the behavior of N for fixed frequency, uj. Since uj = 
ujqc^^i ^"^^^ , tf ~y oo also implies ujq 0. From the 
discussion in Appendix \^ we also know that p ^ oo as 
tJo — > 0. Hence we find 



N{t,L0) 



ujp 



at/{2Rs) 

~7r 



t>tf 



(68) 



This is confirmed by our numerical evaluation of as a 
function of time t > tf for several different values of uj 
(see Fig. [2]). 

Therefore the occupation number at any frequency di- 
verges in the infinite time limit when backreaction is not 
taken into account. This implies that backreaction due 
to particle creation will have important consequences for 
gravitational collapse. 

We have also numerically evaluated the spectrum of 
mode occupation numbers at any finite time and show 
the results in Fig. [3] for several different values of t. The 
similar shapes of the different curves suggest that there 
might be a simple relationship between them. By rescal- 
ing both axes we find that the curves roughly (though 
not completely) collapse into a single curve as shown in 
Fig. 3) Hence, knowing the spectrum at time ti approxi- 
mately gives us the spectrum at all times via 

X-\t) N{t,Lu/X'{t)) = X~\t,) N{t„Q/\'{t,)) (69) 

where, we can determine the function \{t) by considering 
the time variation of N{t,Q), and A' by Eq. The 
result is 



m 

X'{t) 



1 



-t/2Rs 



(70) 
(71) 



We can compare the curve in Fig. 3] with the occupa- 
tion numbers for the Planck distribution 



Np{uj) 



1 



1 



(72) 



where (3 is the inverse temperature. It is clear that the 
spectrum of occupation numbers is non-thermal. In par- 
ticular, there is no singularity in A'^ at w = at finite 
time, there are oscillations in A^, and the rescaling law of 
Eq. ((69)) is not applicable to a thermal distribution. As 
t ^ oo, the peak at w = does diverge and the distri- 
bution becomes more and more thermal. Even at finite 
times, at small frequencies 



1 



Np{uj « /r ^) « — 



(73) 



and the rescaling law amounts to rescaling the tempera- 
ture by a factor AA'. 

Now, from Eq. ([51)) . since the time derivative of the 
wavefunction on the right-hand side is with respect to rj, 
UJ is the mode frequency with respect to rj and not with 
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wRs=l,5, 10,20 




FIG. 2: A'' versus t/Rs for various fixed values of uRs- The 
curves are lower for higher ujRs- At late times the behavior 
is given by Eq. (|^ . 



t/Rs =6,7,8 





t/Rs =6,7,8 



100 



300 



FIG. 4: The same as Fig. |3]but with the axes rescaled as in 
Eq. (|69p . This graph shows that the spectrum at different 
times is approximately self-similar. 

Ln(l + 1/N) t/Rs =8 
0.4 



0.2 



100 



300 
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ojRfi 



FIG. 3: A'' versus ujRs for various fixed values of t/Rs- The 
occupation number at any frequency grows as t/Rs increases. 



respect to time t. Eq. ([55]) tells us that the frequency in 
i is (1 — Rs/R) times the frequency in r], and at time tf, 
this imphes 



UJ 



it) 



(74) 



where the superscript (t) on u; refers to the fact that this 
frequency is with respect to time t. This rcscaling of the 
frequency imphes that the temperature for the asymp- 
totic observer (with time coordinate t) can be obtained 
by find the "best fit temperature" (3~^ and then rescaling 
by (1 — Rs/R)- So the temperature seen by the asymp- 
totic observer is 



T ■ 



~tj/Rs 



(75) 



(The temperature T is not to be confused with the time 
coordinate within the spherical domain wall, also denoted 
by T in Sec. [TTl) By using the scaling in Eq. ([7T|) . it is 
easy to see that /3^^ grows as e+*^/^^^ at late times and 
so T is constant. We can fit a thermal spectrum to the 
collapsed spectrum of Fig. 01 as shown in Fig.Oto obtain 



0.19 



2.4 



Rs iirRs 



= 2-ATh 



(76) 



where Th = l/inRs ~ -08/ Rs is the Hawking tempera- 
ture. Since there is ambiguity in fitting the non-thermal 
spectrum by a thermal distribution, we can only say that 
the constant temperature, T, and the Hawking temper- 
ature are of comparable magnitude. 



FIG. 5: Ln(l + 1/N) versus ljRs for t = 8Rs- The dashed 
line shows Ln(l -I- 1/Np) versus uuRs where Np is a Planck 
distribution. The slope gives f3~^ and the temperature in 
Eq. (ESI). 



The occupation number N{t, u) can be related to the 
asymptotic flux of radiation following standard proce- 
dures (e.g. Chapter 8 of Ref. [l|) and will result in the 
usual greybody factors. 

We thus see that in the context of the Schrodingcr 
formalism there is evidence of Hawking-like, but non- 
thermal radiation emitted during gravitational collapse 
before any event horizon is formed. There are several 
possible sources that one can imagine for this radiation, 
including radiation due to a time-dependent metric, and 
also Hawking emission Q . Since the Schrodinger method 
in principle accounts for all such sources of radiation, it 
is worthwhile reexamining the original Hawking calcula- 
tion, done using the Heisenberg picture and Bogolubov 
machinery, in the context of our above results. 



VI. HAWKING'S CALCULATION 

In Hawking's pioneering paper Q , he considered a col- 
lapsing body. By matching asymptotic field operators, 
he could find the Bogolubov coefficients, and then the 
particle emission rate. The result is the famous Hawking 
thermal radiation at temperature 



Th 



K 

2^ 



(77) 



where k — \/2Rs is the surface gravity. 
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Since Hawking radiation is calculated in the t — > oo 
limit (asymptotic field operators), the result does not 
provide an answer to our original question: what will 
an experimentalist observe at a finite time? So we must 
re-calculate the radiation from a collapsing domain wall 
which is close to, but still larger than, the Schwarzschild 
radius. Stated in a slightly different way - does the exper- 
imentalist sec Hawking radiation before the event horizon 
is formed? 

As Hawking showed, the mode functions of a massless 
scalar field in the black hole spacetime have a "phase 
pile-up" near the event horizon [9|. In other words, if we 
retrace the mode functions from back in time up to 
J'~ 1 the phase of the mode function diverges on J^^ at 
the point vq in Fig. [6l where the coordinate v is defined 

by 




w = i + r + i?5 In 



Rs 



- 1 



(78) 



The radial part of the ingoing mode functions on ,_f 
are (Eq. (2.11) of 0) 



(79) 



The relevant part of the outgoing mode function at 
frequency lo when extended back to J"" is given in 
Eq. (2.18) of [il 



, exp — z — in — 

^/2^r V ^ \ CD 



V < Vq (80) 



and zero for v > vq, where , C, D are constants, and 
K ~ l/2Rs- The expression in Eq. (|80p is only valid for 
small Vq — V, and for large to' (geometrical optics limit). 

(2) ~ 

The overlaps of plj with fi^i and fu,' determine the 
Bogolubov coefficients. This is equivalent to taking the 

(2) 

Fourier transform of Pu ■ Following Hawking's calcu- 
lation, the Bogolubov coefficients for large uj' are (see 
Eq. (2.19), (2.20) of Q; also see H]) 



2tt 



xF ( 1 ) i-iuj') 



/ \ — l+'icj / K 



(2) 



(2) 



(81) 
(82) 



("2") 

Even though the expression for a^^/ is only valid for 
large uj' ^ Hawking argues on analyticity grounds that the 

singularity at w' = should be present. So to obtain 
it becomes necessary to go around the pole at cj' = to 
negative values of lo' . The choice of deformation of the 
contour around the pole is determined on the grounds of 
analyticity, and the result is 



(2) 



cxp 



K / 



„(2) 



(83) 



FIG. 6: Spacetime of a blackhole with null rays originating 
at and going to The last ray that makes it to J^^ 

is emitted at vq. An observer, O, far from the collapsing wall 
will attempt to detect a flux of radiation over a finite but large 
time interval. The last ray to get to the observer originates 
at at Wi < and arrives at ,J^^ a.t u = u\. We are 
interested in finding the particle flux in the section of 
between the points marked and ui. 



' < 



cu' > 



branch cut 

— r^^^ 

Vq - Vi 



FIG. 7: The integration contour in the calculation of the Bo- 
golubov coefficients runs from v = — oo to d = in the com- 
plex 11 plane. In Hawking's calculation vo — vi = 0, and the 
branch cut starts at the origin. For cu' < 0, the integration 
contour is rotated to the upper imaginary axis, and for cj' > 
to the negative imaginary axis. In Hawking's case, this relates 
the Bogolubov coefficients and P^i^i as pedagogically de- 
scribed in Ref. [^. (As pointed out to us by F. Dowker, care 
is required in comparing the calculations in 0] and since 
Hawking considers modes e^'" " while Townsend considers 
g-it^ u -yy^g fQiiQ^ing Hawking's calculation.) In our case, 
these very rotations can also be done. However, the branch 
cut starts at vo — I'l > and the simple relation between 
the Bogolubov coefficients needed for thermal emission is not 
obtained. 
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From here, the calculation of the thermal flux of Hawking 
radiation follows. 

Now consider an observer who only sees the collapsing 
object for a finite time (see Fig. [6]). The last ray detected 
by such an observer emerges from at v = Vi < vq. 
For this observer, the phase of the mode functions have 
a tendency to pile-up but there is no divergence as in 
Eq. ([50)) because v < vi < Vq- As far as this observer 
is concerned, the behavior in Eq. ([50)) holds for v < vi, 
while for v > vi the back-tracked mode functions vanish. 

(2) 

The Fourier transform of pi, now gives the Bogolubov 
coefRcients the following uj' dependence 



(2) 



dv exp 



-iuj V — i — ln(i'o — v) 



(84) 



Following Ref. Q, for w' > we rotate the integration 
contour to the negative imaginary axis (v — s- ~iy) and 
for uj' < to the positive imaginary axis (see Fig. [7]). 
Simple manipulations for > then give 



where 5 = {vq — ui) and = tan ^{S/y). Similarly, for 
w' < we get 



„(2) 



□ (2) 



J2) 



dy e-|"'ly+'^e/«(y2 ^ ^2yu:/2K 



^u}(~iA}') ^^"^ above expressions yield both 



Since (3^ 

Bogolubov coefficients 

The crucial difference between Hawking's asymptotic 
result and the finite time result is the factor exp(±w(?/«;) 
within the integral. Because of this factor, the relation in 
Eq. ()83p does not hold and thermality is lost. However, 
if this extra factor is nearly unity, we can expect the 
spectrum to be nearly thermal. The integral is cut-off 
exponentially for y > l/oj' and hence we estimate that 
the spectrum will be nearly thermal provided ujuj'S/k ^ 
1. Hence the spectrum is thermal at low frequencies and 
gets closer to being thermal as time goes on {d — > 0), 
both of which seem plausible on physical grounds. 

It is difficult however, to go beyond these qualitative 
statements in attempting to compare our results with 
what one might derive in the Hawking approach, in par- 
ticular to determine possibly how much of the effect we 
obtain might be due to Hawking radiation, as opposed 
to particle creation by a changing metric. This is be- 
cause the spectrum depends on a sum over all uj' , while 
the Hawking analysis is done in the geometrical optics 
limit, at large frequencies. Hence to find the spectrum in 
this approach, we need a more complete solution to the 
equations of motion for all the modes of the scalar field 
in the domain wall background. Such solutions are more 
difficult to obtain (as described in [l| for example). 



VII. INFALLING OBSERVER 

So far we have considered the wall collapse from the 
point of view of an asymptotic observer. From the point 
of view of an infalling observer, the time coordinate is r 
of Sec. [IT] and the collapse appears to proceed differently. 
For example, if we ignore radiation, the classical equation 
of motion can be written from the conservation of M in 
Eq. ([T7)) . Then, as the wall approaches the Schwarzschild 
radius. 



M 



- 1 



(85) 



The right-hand side is a non-zero constant, implying that 
the wall is collapsing with constant velocity in the r co- 
ordinate. This shows that the collapse into a black hole 
occurs in a finite time interval for the infalling observer. 
Further, Hawking has argued 0] that the infalling ob- 
server does not detect significant Hawking radiation since 
the emission is dominantly at low frequencies compared 
to l/Rs, while the infalling observer can only have local 
detectors of size less than Rs- Thus the infalling observer 
would appear to see event horizon formation in a finite 
time, with no significant radiation emanating from the 
black hole. 

These paradoxical views of the asymptotic and in- 
falling observers need to be reconciled, and the con- 
ventional way to reconcile them is summarized in the 
spacetime diagram of an evaporating black hole shown 
in Fig. |8| The diagram is drawn so that the asymptotic 
observer sees evaporation in a finite time and the infalling 
observer falls into the black hole in a finite time also. 

We have to note that the diagram in Fig. |5| does not 
follow from a rigorous solution to the problem of radia- 
tion from a collapsing object with backreaction included. 
There are some analyses of this problem in (1 + 1)- 
dimensional models [TO, [HI whose connection with the 
(3 -I- l)-dimensional problem is unclear (e.g. [l2|)- Thus, 
the diagram in Fig. |5| is a conjectured diagram that is 
widely used in the literature. While this diagram may 
well be the correct one once the full problem of gravi- 
tational collapse with backreaction is solved, we have to 
emphasize that it also has some puzzling features that 
indicate that it may not be the best conjecture to make 
in the absence of a backreaction analysis. 

The conventionally drawn spacetime of an evaporating 
black hole has features that are not consistent with our 
findings. Since the asymptotic observer sees Hawking- 
like radiation from the collapsing wall prior to event hori- 
zon formation, the mass of the collapsing wall must be 
decreasing, and at the point denoted by F in Fig. |5| the 
entire energy of the wall has been radiated to J^+. How- 
ever, in the spacetime of Fig. [51 it is at precisely this 
instant that the asymptotic observer sees infalling ob- 
jects disappear into the event horizon, even though there 
is nothing left of the collapsing wall to form the singu- 
larity. A spacetime region such as the triangular region 



11 




i 



FIG. 8: The conventional spacetime diagram for an evaporat- 
ing black hole. The observer A will register a flux of quantum 
radiation even during collapse, and will be able to account for 
the entire energy of the shell by the time he/she gets to the 
line EF. From this point on, A will conclude that there is no 
energy left in the region of the collapsing shell. Yet A will see 
objects and other observers (such as I) disappear into what 
can at most be a Planck scale object and this is a puzzling 
feature of this picture. 



behind the event horizon only seems reasonable if not 
all of the collapsing shell energy has been lost to 
up to the point F, and there is some energy-momentum 
source left behind to crunch up in the singularity. Also, 
if the spacetime near the event horizon is described by 
the Schwarzschild metric, there is infinite gravitational 
redshfit of signals escaping to infinity, while the diagram 
shows that signals escape to infinity in a finite time. Fi- 
nally, as is well known, the diagram in Fig. [S] also gives 
rise to the information loss paradox. While these features 
of the diagram in Fig. [5] arc not inconceivable, they are 
sufficiently strange as to cast doubt on the validity of the 
picture. 

Instead it may happen that the true event horizon 
never forms in a gravitational collapse. We saw that an 
outside observer never sees formation of a horizon in fi- 
nite time, not even in the full quantum treatment. What 
about an infalling observer? As in Hawking's case, the 
infalling observer does not see radiation, but this is due 
to size limitations of his detectors. The mode occupation 
numbers we have calculated will also be the mode occu- 
pation numbers that the infalling observer will calculate, 
even if they be associated with frequency modes that he 
cannot personally detect. The infalling observer never 
crosses an event horizon, not because it takes an infinite 
time, but because there is no event horizon to cross. As 
the infalling observer gets closer to the collapsing wall, 
the wall shrinks due to radiation back-reaction, evaporat- 
ing before an event horizon can form. The evaporation 
appears mysterious to the infalling observer since his de- 
tectors don't register any emission from the collapsing 



wall. Yet he reconciles the absence of radiation with the 
evaporation as being due to a limitation of the frequency 
range of his detectors. Both he and the asymptotic ob- 
server would then agree that the spacetime diagram for 
an evaporating black hole is as shown in Fig. [9l In this 
picture a global event horizon and singularity never form. 
A trapped surface (from within which light cannot es- 
cape) may exist temporarily, but after all of the mass is 
radiated, the trapped surface disappears and light gets 
released to infinity. 

The spacetime picture that we are advocating is similar 
to that described in Refs. [H, [l^ &^<i, more recently, 
Refs. [ll,[ll,[l3. 



VIII. DISCUSSION 

In this paper we have studied the collapse of a 
gravitating spherical domain wall using the functional 
Schrodinger equation. We would like to clearly delineate 
our analysis of the collapse and the emitted non-thermal 
quantum radiation from the interpretational issues about 
the formation of an event horizon. 

First, we studied the collapse of a gravitating spher- 
ical domain in both classical and quantum theory, ig- 
noring any evaporative processes. It has been suggested 
in the literature that quantum fluctuations can cause 
the collapse and formation of a black hole in a finite 
(Schwarzschild) time However, our results show that 
this is not the case and the horizon does not form in a 
finite time even in the full quantum treatment. 

Then we studied radiation from the collapsing shell as 
seen by the asymptotic observer. In the process of gravi- 
tational collapse, there are two, perhaps related, sources 
of radiation: first is the radiation from particle creation 
in the changing gravitational field of the collapsing ball, 
and second may be Hawking-like radiation due to a mis- 
match of vacua at early and late times. The functional 
Schrodinger analysis takes all such sources into account 
and therefore gives the total particle production. We 
have found a non-thermal distribution of particle occu- 
pation numbers, with departures from thcrmality as illus- 
trated in Fig. [3] and discussed toward the end of Sec. |Vl 
In a limited range of frequencies, the spectrum is approxi- 
mately thermal and the temperature fitted in a restricted 
range of frequencies is constant and roughly equal to the 
Hawking temperature l/AnRs- The radiation becomes 
thermal in the entire range of frequencies only in the 
limit t —f oo, i.e. when the horizon is formed. Further, 
the mode occupation number diverges in the infinite time 
limit; if the backreaction is neglected {i.e. the back- 
ground is held fixed). Since an outside observer never 
sees formation of a horizon in a finite time, radiation 
observed by him is never quite thermal. (Non-thermal 
features also get greatly amplified once the background 
is also treated quantum mechanically [3.) This non- 
thermal radiation has strong implications for the infor- 
mation loss paradox since it can carry information about 
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the collapsing matter. 

Without a rigorous calculation that includes backrc- 
action, one can not give a definite answer to the final 
fate of a collapsing object. It may happen that the dia- 
gram in Fig. [5] is correct and some radical and elaborate 
solutions to the problems we mentioned in Sec. IVIII are 
needed. However, one can imagine an alternative pic- 
ture, different from the one in Fig.[8l which seems to have 
fewer problems, and that is that an event horizon never 
forms. Since the mass of the shell is decreasing during the 
collapse, the shell will be chasing its own Schwarzschild 
radius, and the question is whether the shell will catch up 
to its own Schwarzschild radius or completely evaporate 
before that happens [T3 |. 

With backreaction included, the radiation should lead 
to a continual reduction of the Schwarzschild radius, Rs, 
occurring in the Ipser-Sikivie metric (see Sec. |TT|. Then, 
as seen by the asymptotic observer, one of two possi- 
bilities occurs: either the collapsing domain wall evapo- 
rates and Rs ^ in a finite time, or else backreaction 
causes the radiation rate to slow down and vanish in a 
finite time. This latter possibility is unlikely, as our esti- 
mates suggest that the rate of emission increases as Rs 
decreases [23|. We therefore conjecture that the backre- 
action due to particle production will cause the collapsing 
domain wall spacetime to completely evaporate in a finite 
time. In this case, the spacetime can either be as given 
in Fig. [51 or have the same global spacetime structure as 
Minkowski space, as shown in Fig. [5] If the latter pic- 
ture is correct, it also means that the infalling observer 
will not encounter an event horizon, because this feature 
is simply absent from the spacetime. Another way to 
see this is to note that the causal relation between two 
events is the same for all observers. Hence if the asymp- 
totic observer sees a signal from an infalling object after 
he sees the last radiation ray emitted by the evaporating 
wall, this will also be the sequence of signals seen by the 
infalling observer. As discussed in Sec. IVIIi the infalling 
observer would expect to see an intense burst of radiation 
as the wall approaches the Schwarzchild radius, but can 
fail to do so because his detectors are too small to detect 
the emitted range of frequencies. 

In the absence of an exact backreaction calculation, 
we also have to allow for the possibility that a value of 
the critical mass exists above which Fig. [S] applies and 
below which Fig. [5] holds. Also, as discussed by Hawk- 
ing [l^, the question of "whether a black hole forms" 
is not sharp enough and may not make sense in the full 
quantum theory since all of the measurements are made 
by an asymptotic observer at infinity, while a collapsing 
object exists for a finite time and disappears by emit- 
ting radiation in the strong field region in the middle. 
An asymptotic observer can never be sure if a black hole 
formed because of underlying quantum uncertainty [T^ . 

The broad picture we have obtained is consistent with 
that proposed in Refs. [13, [l3|j though there are differ- 
ences in the analysis and the conclusions. In particular, 
we find a non-thermal spectrum whereas Gerlach argues 
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FIG. 9: The spacetime of a collapsing domain wall. During 
collapse the wall emits non-thermal (quasi-Hawking) radia- 
tion as depicted by the arrows. Our calculations indicate 
that the total energy flux between the point i° to some point 
indicated by F is equal to the energy of the initial domain 
wall. Hence we conjecture that the domain wall evaporates 
completely at point D. Between F and i^, there is no radia- 
tion flux arriving at The event horizon and singularity 
present in the customary treatment are not formed and the 
spacetime structure is the same as that of Minkowski space- 
time. 



for thermality. Our picture also supports the interpre- 
tation of Hawking radiation given in Ref. Q whereby 
particles are created during the process of gravitational 
collapse and are then radiated slowly to form what we 
call Hawking radiation. We have indeed found particle 
production during the collapse but the radiation is not 
quite thermal. It is only in the frequency range where 
the occupation number spectrum can be approximated 
by T/uj (Eq. ([75)1 ) that thermality holds at finite time. 
Also note that the non-thermality wc find is in the mode 
occupation numbers. Propagation of the radiation in the 
background metric will cause further non-thermality due 
to greybody factors. 

If we live in a world of low scale gravity, the colli- 
sion of particles in high energy accelerators will lead to 
a situation where the particles are in a continual state of 
gravitational collapse from which non-thermal radiation 
is being emitted. The life-time of such a state can be es- 
timated once we know the details of the radiation more 
precisely from an analysis which includes backreaction. 
However, on dimensional grounds, Hawking's estimate 
for the lifetime of a black hole (~ i?|/G) may well apply 
to the colliding particles as well. 

In reality the collapse is further complicated by the 
fact that the collapsing object is not kept in isolation 
and there are external forces that can disrupt the collapse 
at any point in time. From the perspective of potential 
information loss, note that any infalling encyclopedias 
can be returned to the asymptotic observer if the collapse 
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is disrupted at any time, as it could be, for example by 
a bomb set to go off at some late, but finite time. Most 
importantly, since we calculate that the radiation emitted 
during the gravitational collapse is never truly thermal, 
the classic information loss issue in black holes should, in 
this case, be a non-problem for the asymptotic observer 

Our primary result, that no event horizon forms in 
gravitational collapse as seen by an asymptotic observer 
is suggestive of the possibility of using the number of 
local event horizons to classify and divide Hilbcrt space 
into superselection sectors, labeled by the number of lo- 
cal event horizons. Our result suggests that no opera- 
tor could increase the number of event horizons, but the 
possibility of reducing the number of pre-existing primor- 
dial event horizons is not so clear and would require that 
Hawking radiation not cause any primordial black hole 
event horizons to evaporate completely. 

Our conclusions have been derived on the basis of a 
number of assumptions which we now discuss. The first 
is the truncation of supcrspacc to minisupcrspacc. We 
have only included spherically symmetric field configu- 
rations. Even then, the metric is restricted to be the 
classical solution sourccd by a spherical domain wall. A 
more general analysis would include more metric degrees 
of freedom, though it is hard to see how this would make a 
difference to our conclusion. Similarly, we have restricted 
ourselves to a zero thickness domain wall. A more gen- 
eral analysis would allow for a thick wall. Finally, the 
Wheclcr-dc Witt formalism as we have used it, does not 
allow for the creation and annihilation of domain walls 
("third quantization"). Perhaps third quantization could 
allow for the spontaneous creation of a black hole and the 
annihilation of the wall, effectively leading to black hole 
formation. A fourth possibility is that our Lagrangian 
breaks down near the Schwarzschild horizon and "quan- 
tum gravity" effects become important. This is usually 
thought not to be the case since the spacetime curvature 
near the horizon is small for large black holes. 

Perhaps the most serious drawback of our analysis is 
that it does not include backreaction on the gravita- 
tional collapse due to radiation. While we do not expect 
such inclusion to alter our conclusions regarding the non- 
existence of event horizons for asymptotic observers, we 
are currently exploring ways to extend our treatments to 
include backreaction. 

No theoretical idea is complete without the possibility 
of experimental verification and so it is important to ask 
if the picture we have developed in this paper can also 
be tested experimentally. We have already mentioned 
the relevance of our conclusions to black hole produc- 
tion in particle accelerators provided low scale gravity is 
correct. However, there is an even more accessible ex- 
perimental system where these theoretical ideas can be 
put to the test. These are condensed matter systems in 
which sonic black holes (dumbholes) may exist [2^. It 
is very hard to realize a dumbholc in the laboratory for 
various experimental reasons; the closest known realiza- 



tion seems to be the propagating Hc-3 AB interface in 
the experiment of Ref . [2l| as discussed in . Yet the 
crucial aspect of our work in this paper is that there is 
no need to produce a dumhhole in order to see acous- 
tic "pre- Hawking" radiation. The process of collapse to- 
ward a dumbhole will give off radiation. This is also the 
conclusion of Ref. [l^l though the details of the analy- 
sis and conclusions are different - for example, we find 
non-thermal emission whereas these authors claim ther- 
mal emission with a modified temperature that is lower 
than the Hawking temperature. In any case, it should be 
much easier to do experiments in the laboratory that do 
not go all the way to forming a dumbhole, and this could 
be an ideal arena to test pre-Hawking radiation. 

Our conclusions are important not only for the gen- 
eral issue of the breakdown of unitarity via information 
loss, but also for more general studies of black hole for- 
mation, whether they be in the context of astrophysics 
{e.g. galactic black holes) or in future accelerator ex- 
periments. In all these situations, we are asymptotic 
observers watching the gravitational collapse of matter, 
and we may never see effects associated with a black hole 
event horizon. Only effects occurring during the gravita- 
tional collapse itself appear to be visible. 
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APPENDIX A: p EQUATION 

In the range t < 0, w is a constant and the solution to 
Eq. jeO]) is 

Piv) = (Al) 

In the range < t < tf, we do not have an analytic 
solution but we can derive certain useful properties. First 
note that in terms of rj 
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Then the equation for p after rescahngs can be written 
as: 



s) 



f 



1 



1 - 7/ /3 



(A3) 



where rj' = rj/Rs, f — y^uJop- The boundary conditions 
are 



/(0).l, ^M.o 



(A4) 



The last term with the becomes singular as / ^ 0. 
Let us consider another equation with the 1/ replaced 
by something better behaved. For example, 



dif 



1 -r/ 



9 



with boundary conditions 



(A5) 



(A6) 



Eq. (jAsp implies that is monotonically decreasing 

as long as g{'q') > 0. Furthermore, it is decreasing faster 
than the solution for / as long as / < 1, since the 
term in Eq. (jA3|) is a larger "repulsive" force than the g 
term in the Eq. (KE\i . So 



9iv') < fiv') 



(A7) 



for all rj' such that 5(77') > 0. 

Eqs. (jASP with initial conditions (|A6[) can be solved 
in terms of degenerate hypergeometric functions. For us, 
the important point is that the solution for g is positive 
for all ij' and, in particular, ^(l) > for all the values of 
ujqRs that we have checked. Therefore is positive, 

at least for a wide range of luqRs- 

Let /i = /(I) 7^ 0. Then the equation for / can be 
expanded near ?/' = !. 



dif 

This shows that 
drj' ^ 



-{ujqRs) 



h 



1 - '/ ft 



{uooRsY fi ln(l - 77') -00 



(A8) 



(A9) 



as 77' "> 1. 

Hence ^(77 = Rs) is strictly positive and finite while 
Prii'T] = Rs) = —00 for finite and non-zero loq. Since 
/ = y/^p, and / — > 1 for loq ^ 0, we also see that 
p ^ 00 and Pjj a.s luq ^ 0. 

In the range tf < lo is a, constant. However, the 
solution for p is not a constant, unlike in the range t < 0, 
since the constant solution 1/ yjLo(tf) does not necessarily 
match up with p(tf—) to ensure a continuous solution. 
Yet it is easy to check that in this region = and 
so there is no change in the occupation numbers. So we 
need only find N{tf — ^uj) to determine N{t —f 00, (D). 



APPENDIX B: NUMBER OF PARTICLES 
RADIATED AS A FUNCTION OF TIME 

We use the simple harmonic oscillator basis states but 
at a frequency w to keep track of the different w's in the 
calculation. To evaluate the occupation numbers at time 
t > tf, we need only set lu = uj(tf). So 



HniV^b) (Bl) 



where Hn are Hermite polynomials. Then Eq. ([65]) to- 
gether with Eq. ([59)) gives 



where 



1 ^ 


r 




TT^Ujp-^ ^ 




V2"n! J 


1 ^ 


r 


ia 

r In 


TT^Ojp'^ ^ 




V2"n! 


P = 


1 - 


Q\p 



"n! J 



(B2) 



(B3) 



To find /„ consider the corresponding integral over the 
generating function for the Hermite polynomials 



J{z) 



Since 



z=0 



Therefore 

Since 

Hn{0) = 



(B4) 



(B5) 



(B6) 



(B7) 



ji = even (B8) 



and Hn{0) = for odd n, we find the coefficients c„ for 
even values of n, 



(-l)"/2e*" / 2 



((1-^2)1/4 V P 

For odd n, c„ = 0. 



( 1 - ^ I ^ J (B9) 
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Next wc find the number of particles produced. Let 

2 



X 



1 



P 



(BIO) 



Tlren 



n— even 

2 



± ^ (n-1)!! „ 



I" I n— even 

2 d 1 



2 X^ 

^\P\ (l-x2)3/2 



(Bll) 



Inserting the expressions for x a-nd P, leads to 



N{t,oj) 



up 

7f 



Prj_ 

ujp 



(B12) 



In summary, we have found the occupation number 
of modes as a function of p which is a function of time 
as given by the non- linear differential equation Eq. (j60p . 
The equation connecting p and time t has only been 
solved numerically but we have discussed the behavior 
of p and p,) as ry — > Rs [t oo) in Appendix [Al 
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